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Abstract 

We study the neutron electric dipole moment in the presence of the CP-violating 
operators up to the dimension five in terms of the QCD sum rules. It is found that 
the OPE calculation is robust when exploiting a particular interpolating field for 
neutron, while there exist some uncertainties on the phenomenological side. By us- 
ing input parameters obtained from the lattice calculation, we derive a conservative 
limit for the contributions of the CP violating operators. We also show the detail 
of the derivation of the sum rules. 



1 Introduction 



A variety of experimental efforts [T] has precisely determined the elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix [2], |3], which is a source of CP violation in the 
Standard Model (SM). All of CP- violating processes observed ever are well-explained 
in terms of the single physical phase in the CKM matrix. The SM, which is based on the 
SU{3)c X SU{2)l X U{1)y gauge symmetry, allows another CP- violating interaction: the 
9 term in the QCD sector. The CP-violating phenomena caused by the interaction are, 
however, quite different from those induced by the CKM phase; the QCD 9 term gives 
rise to the CP violation in the flavor- conserving processes, while the CKM phase induces 
the CP violation in the flavor-changing ones. Furthermore, TeV-scale physics beyond the 
SM, such as the Minimal Supersymmetric Standard Model (MSSM), often provides other 
CP-violating sources. In fact, additional CP- violating interactions are necessary from the 
cosmological point of view, since the observed baryon asymmetry in the Universe is not 
explained within the SM interactions. 

The electric dipole moment (EDM) of neutron is one of the physical quantities that 
are quite sensitive to the CP violation in the flavor-conserving interaction. Since there 
has been no experimental evidence for its existence so far, a severe constraint is imposed 
on the CP-violating interactions. The currently most stringent limit for the neutron EDM 
is given by the Institut Laue-Langevin (ILL) experiment [4j: 

\dn\ < 2.9 X 10"^^ e cm (90% C.L.) . (1) 

Moreover, several experimental projects which use ultra cold neutrons are now under 
development and expected to have much improved sensitivities. For example, the nEDM 
collaboration at the Paul Scherrer Institute (PSI) [5] plans to deliver a sensitivity of ~ 
5 X 10~^^ e cm, and eventually to reach into the regime of 10~^^ e cm. Similar sensitivities 
are expected to be achieved by the nEDM Collaboration at the Spallation Neutron Source 
(SNS) at the U.S. Oak Ridge National Laboratory [6j, the CryoEDM experiment [7], 
the NOP Collaboration at J-PARC [8J, and the experiment at KEK-RCNP-TRIUMF 
[H]. Such high sensitivities provide an opportunity to probe the flavor-conserving CP- 
violating interactions in the TeV-scale physics beyond the SM. Furthermore, we may 
probe the flavor violation in the new physic indirectly. Even if the new flavor- violating 
interactions are introduced in the new physics, the relative CP phase between them and 
the CKM matrix may contribute to the EDM [TU] . 

In order to translate the experimental limits for the neutron EDM into constraints on 
the CP violation on the Lagrangian at parton level, one needs to obtain a relation between 
these two quantities. There are some attempts to derive the relation based on the naive 
dimensional analysis, the chiral perturbation theory, and the QCD sum rules, though they 
are considered to have large uncertainties. It is ultimately desired that the lattice QCD 
simulation would evaluate it in future. There has been discussion of evaluation of neutron 
EDM with lattice simulation [TT] . 

In this work, we evaluate the neutron EDM with the QCD sum rules [12], including the 
CP-violating operators up to the dimension five. It is considered that the QCD sum rules 
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allow us to derive the relation more systematically than the naive dimensional analysis 
and the chiral perturbation theory [13]. Similar attempts have been already made in 
the previous works, e.g., in a series of papers by M. Pospelov and A. Ritz [HI [15] and 
references therein. We also derive the sum rules for the neutron EDM, while we use the 
lattice QCD simulation result for the low-energy constant in the numerical evaluation of 
the neutron EDM. It is found that this gives more conservative estimate than carrying 
out all of the evaluation within the framework of the QCD sum rules. This approach 
provides a way of eliminating theoretical errors from the calculation, while there still 
remains uncertainty resulting from the QCD sum rule technique itself. 

This paper is organized as follows. In Sec. [2l we review the CP- violating interactions at 
parton level up to the dimension five. From Sec. [21 the analysis of the neutron EDM with 
the QCD sum rules starts. In Sec. [3]we discuss phenomenological aspects of the correlator 
of the interpolating field to neutron, and, in Sec. [U show the properties of the neutron- 
interpolating field. In Sec. [5l the quark propagators are derived on the CP-violating and 
electromagnetic background. They are used to evaluate the operator product expansion 
(OPE) for the correlator in Sec. [6] The sum rules for the neutron EDM are derived in 
Sec. [3 We found that there is a difference between results in Refs. [TU [15] and ours. 
In Sec. [HI we extract the low-energy constant from the lattice QCD simulation result. In 
Sec. [9] our numerical results for the neutron EDM are derived. In Sec. [10] the neutron 
EDM is discussed assuming the Peccei-Quinn symmetry solves the strong CP problem 
[16] . Section [TT] is devoted to conclusion and discussion. 

In Appendix, we show some useful formulae to derive the quark condensates in the 
CP-violating background. In Appendix \^ we estimate the effect of the CP-violating 
interactions on the generic quark bi-linear condensate (0|grg|0), with P a 4 x 4 constant 
matrix, as well as on the quark and gluon background fields. In Appendix [B] validity of 
usage of the classical equations of motion of quarks in evaluation of the quark condensates 
is discussed. In Appendix [O the Wilson-line operators for the quark fields are discussed 
in the Fock-Schwinger gauge. 

2 Effective Lagrangian 

Let us first express the flavor-conserving CP-violating terms in the low-energy effective 
Lagrangian for the system consisting of light quarks and gluon. We include all of the 
CP-violating operators up to the dimension five: 

^ — ' ovr ^ 

q=u,d,s 

- ^ X] ^^^'^^ ■ - ^ 5Z ^<i^9s{G ■ a)-i^q . (2) 

q=u,d,s q=u,d,s 

Here, rrig represents the quark masses, and G^^ are the electromagnetic and gluon field 
strength tensors, gs is the strong coupling constant (a^ = gl/An), F ■ a = F^^cr'^'^, G • a = 
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G^^a'^'T^, and G^^ = le^^^p^G^P" with e^i^s = +i. denotes the generators in the 
SU{3)c algebra. The second, third and forth terms in Eq. are called the effective QCD 
6 term, the electric and chromoelectric dipole moments (CEDMs) for quarks, respectively. 
The EDMs and CEDMs for quarks are dimension-five operators, and they are sensitive 
to the TeV-scale physics beyond the SM. The coefficients of the CP-violating operators, 
9g, 9c, dq, and dg, are all assumed to be quite small, and we keep only the terms up to 
the first order of these parameters. 

The first two terms in Eq. (|2]) are mutually related by the chiral rotation. Consider 
the following infinitesimal chiral rotation: 

q^q' = (I- iepq-i^) q , (3) 

where e is an infinitesimal real constant and Pq are certain parameters for each quark. 
The Noether current associated with the transformation is given as 

= 5Z /'9^7A.75g • (4) 

q=u,d,s 

The divergence of this current does not vanish. Instead, 

- X] P'Mi^^^ + dqq9s{G ■a)q] . (5) 

Hereafter we choose pg as 

Pq = Oq/OQ, Oq^ 5^ 0, . (6) 

q=u,d,s 

Then, if we take the infinitesimal parameter in Eq. ([3]) as e = Oq/2, the Lagrangian in 
Eq. ([2]) varies by 




= ^m,g^^,75g + ^Q-^G;:,G«^% (7) 

1 

which implies that 

q=u,d,s q=u,d,s 

where 9 = 9g + Oq- 

Therefore, it is found that the 75-mass terms are always reduced to the ordinary ones, 
and it is 9 that is regarded as a physical parameter. Of course, one may in turn rotate 
out the 9 term into the imaginary mass term through an appropriate chiral rotation. 
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In addition, there remains still some arbitrariness in the quark mass phases 6'^, since 
they are redefined into another through an SU{3) chiral rotation. In this article, we 
choose an appropriate set of 9q so that the choice significantly reduces the CP-violating 
contribution to the vacuum expectation values (VEVs) of quark bi-linear. We take the 
condition in Ref. flT] to determine 6q, that is, after the 6 term rotated into the 75-mass 
term, the following relation should be satisfied: 



{n^\C^\M^) = , iM^ = n, K, rj) . 



(9) 



The above condition is evaluated by using the partially conserved axial-vector current 
(PCAC) relations. In the current case, it is sufficient to examine the conditions for 7r° 
and ?7°. By using the PCAC relations, one may readily deduce the conditions for the 
CP- violating parameters from Eq. (|9]): 



-ml{du + dd- 24) . 



(10) 



0{mupu - mdPd) = -ml{du - dd 

1 

OiniuPu + rUdPd - 2msPs) 

In the calculation we parametrize the condensate {qgs{G ■ a)q) as [H] 

{qgs{G ■ a)q) = -mliqq) ■ (11) 
With the relation ^ pg = 1, we then determine the quark mass phases as follows: 



Pu = 

rriu 

Pd = — 
rrid 

Ps = — 



1 + ^ 
26 



rrin 



1 + ^ 
26 



rrin 



1 + 
26 





+ 


du dg 






rrid 




ITLs 


dd — du 


+ 


dd — ds 




ruu 




rus 


dg c/y 


+ 


ds — dd 



where 



mumdrris 



rrid 



ruumd + mdrus + ruums 



(12) 
(13) 



3 Phenomenological behavior of correlator 

The QCD sum rules are based on an analysis of the correlator of interpolating field^]. 
In the method, OPE allows one to consistently separate the long- and short- distance 
contributions to the correlator, and the long-distance contributions are evaluated by con- 
densations of quarks and gluon. By comparing the evaluated correlator with the phe- 
nomenological model, the properties for the low-lying parts of the hadronic spectrum are 

^ There are many review articles about the QCD sum rules. For example, see Refs. [19]. 
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derived. The Borel transformation is applied to the correlator there. In this section, we 
first discuss the phenomenological model for the correlator. 

In the present case, the interpolating field must have the same quantum numbers 
as those of neutron, and it is denoted by rin{x) hereafter. On a background with CP- 
violating sources, the matrix element of the interpolating field between the vacuum and 
the one-particle neutron state is given as 

{Qc/p\Vn{x)\Nc/p{p, s)) = zic/p ■ Unjp{p, s) e-'P-^ , (14) 

where \^cIp) and \NcJp{p, s)) indicate the vacuum and the one-particle neutron state 
on the CP- violating background, respectively. The spinor Un,c/p{p, s) is on-shell neutron 
wave function which satisfies the Dirac equation: 

- m^^c/p ■ e-^'^"^«)M„,c/p(p, s) = 0. (15) 

Here we include a phase factor e"*""'''^ into the mass term, which in general might appear 

1 

as CP is broken in the vacuum. Since Zn,c/p and mn,cjp are both even in terms of the 
CP- violating parameters [HJ, up to the first order of them, 

1 1. 

Z^^cjp = Zi, mn,c/p = rUn , (16) 

where m„ is the mass of neutron and A„ = Zn is the coupling between the physical 
neutron state and the interpolating field without CP-violating sources. Then the solution 
of Eq. (fT5|) turns out to be 

Wn,c/p(p,s) = et-"^«M„(p,s) , (17) 

with UniPi s) an ordinary spinor wave function which satisfies (|^ — mn)Un{p, s) = 0. As a 
result Eq. ( IT^ leads to 

{Qc/p\Vn{x)\Nc^pip, s)) = A„e^""^^ M„(p, s) e"^^'^ . (18) 

The low-energy constant A„ is to be determined later. 

Now we analyze the correlator of the interpolating fields from the phenomenological 
viewpoint. It is defined as 

n(g) = tjd'x e'"-^ {Qc/p\T{Vn{x)VnmWp)F , (19) 

where the subscript F implies that the correlator is evaluated on an electromagnetic field 
background. Our goal is to extract the EDM of neutron from the correlator. The phase 
factor in Eq. f lTH]l . however, causes mixture between electric and magnetic dipole moment 
structures and makes it difficult to pick out only the EDM from the QCD sum rules. So 
we first examine the Lorentz structures of the correlator and select a term independent 
of the phase i.e., chiral invariant. As discussed in Ref. [H], up to the leading order 
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on the background electromagnetic field, the correlator n(g) is estimated by inserting an 
effective vertex such as 

Cn = -\dnN{F ■ a)75iV = ^NF ■ aN . (20) 

Here, A^(= N{x)) denotes the renormalized neutron field which is approximately equiva- 
lent to A~^e~*""'^^''^?7„(a;), and dn is the EDM of neutron. A similar procedure to those in 
Ref. [H] shows that terms with an odd number of Dirac matrices are independent of the 
phase factor q;„, and furthermore, those proportional to {F ■ a, ^} are the unique choice 
in this case. Therefore we only focus on such terms in the following calculation. Then, 
the phenomenological expression of the correlator is found to b^. 

n(p^'^-\q) = ^fiq^){F.aJ} + ... , (21) 
where dots indicate terms with other Lorentz structures and 

with v4(g^) and B{q^) functions which have no pole at = m^. As noted in Ref. |14] . 
since we are effectively dealing with a three-point function, it might be inconsistent to 
parametrize the continuum contribution in terms of a usual ansatz for the spectral function 
with a certain threshold in the QCD sum rules. We just neglect the contribution with 
expecting its significance to be small enough. Furthermore, we assume that the function 
A{q^) has little dependence on g^, and regard it as a constant when we conduct the Borel 
transformation. 



4 Neutron-interpolating field 

In this section we give discussion on choice of the neutron-interpolating field which we use 
for the QCD sum rule calculation. The field must have the same quantum numbers as 
neutron. The most general interpolator for neutron on the ordinary CP-even background 
is parametrized as 

Vni^x) = ji{x) + I3j2ix) , (23) 

where 

ji(x) = 2eabc {d'^{x)C-f5Ub{x)) d^{x) , (24) 

and 

jiix) = 2tahc {dl{x)Cubix)) 'j^dcix) . (25) 

Here the subscripts, a, b, c, denote the color indices and C is the charge conjugation matrix. 
The interpolator ji{x) is often used in lattice simulations. While j2{x) vanishes in the 

^ In the published versions of Refs. [UJ [TS] the coefficient of the double pole in Eq. is different 
from ours by a factor of two, while that in the revised arXiv versions are consistent with ours. 
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non-relativistic limit, it should be included to the whole interpolating field since we deal 
with light quarks. The unphysical parameter /3 is to be fixed later so that the calculation 
is transparent. 

When the calculation is carried out on the CP-violating background, however, the 
interpolating fields include additional components. This point is easily understood when 
one considers the chiral rotation discussed in Sec. [2J As we have seen in Sec. [21 the chiral 
rotation transforms the Lagrangian C into another. The same transformation, in turn, 
changes the interpolators ji{x) and j2{x) into other ones as 

ji{x) - ze[(p„ + pd)ii{x) + Pdi2{x)] , 

j2{x) J2ix) - ie[{pu + pd)i2ix) + pdii{x)] , (26) 

where 

ii{x) = 2eabc{d'^{x)Cub{x))dc{x) , 

i2{x) = 2eabc{d^{x)C-f5Ub{x))j5dc{x) . (27) 

Therefore, with generic CP- violating terms as in Eq. ([2]), the interpolators acquire mixing 
terms as 

jiix) = ji{x) + idi{x) + i5i2{x) , 
i2{x) = j2{x) + iei2{x) + i5ii{x) , 
= ii{x) + ieji{x) + i5j2{x) , 
i2{x) = i2{x) + iej2{x) + i6ji{x) , (28) 

with e and S the small constants which are suppressed by the CP-violating parameters. 
Furthermore, the above expressions are rewritten as 

ji{x) = (1 + i575)[ji(x) +ieii(x)] , 
j2(x) = (l + i6-f5)[j2ix) + iei2{x)] , 
ii{x) = (1 + i5'j5)[ii{x) + ieji{x)] , 

i2{x) = {1 + i6^5)[i2{x) + iej2{x)] , (29) 

because 

nix) = 75j2(a;) , 

i2ix) = 75ji(x) . (30) 

Now that we concentrate on the chiral-invariant structure in the correlator of the neutron- 
interpolating field as discussed in Sec. [3|, the infinitesimal chiral rotation factor (1 + iSj^) 
is ignorable. After all, the neutron-interpolating field which we deal with has a following 
structure: 

Vn{x) = ji{x) + /3j2ix) + ie[ii{x) + /3i2ix)] . (31) 
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The small constant e is determined by the condition that the interpolating field rjnix) 
has a vanishing correlator with the current ^n(x) defined as follows: 

^n{x) = ii{x) + Pi2{x) + ie[j,{x) + Pj2{x)] . (32) 

In what follows, however, we sweep away the contribution of the mixture terms in the 
interpolating field by choosing an appropriate value for the parameter (3. In the subsequent 
sections, we calculate the correlator of ?7„(x) by using the OPE method. The correlator 
is expressed by the sum of the correlators for each component interpolator as 

{ncJp\T{r]^{x)UO)}WF = (Jl,jl)+/3[(jlj2) + (j2jl)]+/3'(j2j2) 

+ie[l5{j2,Ji) + {31,32)15] 

+^e/3[{(jl,Jl),75} + {(j2,j2),75}] 

+2e/32[75(ji,j2) + (j2ji)75] , (33) 

where 

(a, 6) = {n^\T{aix)biO)}\n^)F . (34) 

As discussed in Sec. [31 we focus on parts of the correlators which have the Lorentz struc- 
tures with an odd number of gamma matrices. Such terms anti-commute with 75. Thus, 
in this case, the above expression leads to 

{n^\T{7]4x)T]nm\^^)FU odd = (jljl) +/3[(jlj2) + (j2, Jl)] + (j2, J2) 

+te{l-/3')[{3u32)-{32,3i)H- (35) 

This equation shows that the mixing terms in the interpolating field do not affect the 
correlator if one sets P to be ±1. As we will see later, for our calculation, /3 = +1 is 
an appropriate choice since this choice eliminates the sub-leading terms with infrared 
logarithm, which yield ambiguity due to the infrared cutoflEl. With this choice one may 
simultaneously exclude the contribution of the mixing terms. Thus, we will not calculate 
such mixing contributions with keeping in mind that we will finally take P = +1 when 
we derive the QCD sum rule^. That is to say, we deal with the correlator 

{ncJp\T{r]^{x)Um^<Jp)F = (jiJi) +/3[Oi,j2) + 02,Ji)] +/3'(J2,J2) , (37) 

and after the computation, we set /3 = 

The choice of f3 for the sum rules including only the QCD 9 term is discussed in Ref. [3D]. They 
argue that optimal choice of/3is— 1</3<0 rather than 1, which is consistent with the conventional 
choices favored from a viewpoint of evaluation of A„. The discussion is, however, not applicable to the 
present case since our sum rules contain several unknown parameters. 
^ The neutron- interpolating field for /? = +1 is simply expressed as 

Vn{x) = le'''''{dlC<j,,,dt)<J^''j5Uc ■ (36) 
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5 Quark propagators on the CP- violating background 

When evaluating the correlator flT9|) in the OPE, we need to obtain the quark propagators 
on the CP-violating background with an electromagnetic background field F. They are 
defined as follows: 

[Sl{x)U ^ {^dp\T [qaa{x)q,pm \ndp)F , (38) 
where a and /3 denote spinor indices. Expanding the propagators as 

Kb(a;)]a/3 = ^Ti.^) „ + Xla{^)xlp{^) + [Slb{x)] li photon + [Sli{x)\^p |i gi^on + • • • , 

(39) 

we evaluate each term in x-space. The first term is the free propagator, and the sec- 
ond term describes the correlator of the quark background fields, with x'iai^) ^ classical 
Grassmann field which indicates the quark background field. The third and fourth terms 
represent the propagators including one photon and gluon, respectively. In the derivation 
of the quark propagators we use the classical equations of motion for quark fields given 
as 

ilpq = mq{l + ieq-f5)q + ^ ' ^^^^g + ^ ^i9s{G ■ a)^5q , 

q=u,d,s q=u,d,s 

-iq% = mqq{l + iOq^^) + ^^^(^ " o^)75 + ^g9sq{G ■ (7)75 , (40) 

q=u,d,s q=u,d,s 

where = — iCgA^ — ig^G^T^ and qlp = d^q'j^ + ie^q^ + igsqCfi^T^. The elec- 
tromagnetic and gluon fields are denoted as and G^, respectively, with the quark 
charge. 

The first term in Eq. fl39|) . [S'^ab\x)\ai3-, is readily evaluated by using the equations of 
motion without electromagnetic and gluon background fields. The result is 

gqio), . _ _ _ , . 

\b - (3.2)2 4^23,2 , 

where we keep the terms up to the first order of the quark mass rriq. 

Next, we evaluate the third and fourth terms in Eq. ( 139|) . These terms again are 
obtained from the equations of motion ( 140|) . In this calculation, it is convenient to use 
the Fock-Schwinger gauge [21] for both the electromagnetic and the gluon fields: 

x^A^(x) = x^'G^ix) = . (42) 

In this gauge, the fields are expanded by their field strength, such as 

M^) = + ^x'^x'^iD^F^.m + -L.x-x^:,'^^D^Df,F,,iO)) + ■ ■ ■ . (43) 
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By using the expression, the gauge covariant form of the propagators is obtained as follows: 



•^06(^)11 photon 



te 



327r2 

dg 



1 \ 

Oab 



ab 



X 



imq[l 



z^,75)F-alog(-A?^x2) 



r.2\2 



+ 



ma 



(44) 



SI. 



327r2 
87r2 



— {^,Ga6 ■ 0"} -imq{l -i6q-f5)Gab ■ crlog(-Aj^, 



2 2 



(a;2)2 + 2x2 



(45) 



with = G^^'^T^. Here, a certain infrared (IR) cutoff A//j is introduced in logarithmic 
terms. It was replaced to the quark masses when deriving the propagator from the 
equation motions. However, the contribution to the OPE with small quark momenta 
around the quark masses should not be included so that the IR cutoff is introduced. 

Finally, we translate the quark and gluon background fields into their condensates. 
Here we just give resultant expressions for the relations between them. The details of the 
derivation is presented in Appendix |X1 

A single quark line, Xaa(3^)Xb/3(0), is related with the quark condensate as 



and it is to be expressed in terms of (qq) as follows: 

Sab 



(46) 



iX^,(0) 



^2 + ^^GPgTs)^^ (qq) + -^Sab^a^rrigiqq) 



96 



OnigPgCqX + dq + {K- -i)eqdq 



{F ■ (y4}ap{qq) 



+ -^rriqCqX^abiF ■ a,:/:}^p{qq) 



24 



~ TT^^^gXSab ( F ■ (T75 [1 + ipg^GTs] ) (qq) ■ 



(47) 



Here, Xy ^ind ^ are the parameters for the quark condensates defined as p2] 

{q(^i,uq)F = eqxF^u{qq), (48) 
gs{qG^.T^q)p = eq^F^,{qq) , (49) 



'^9s{qibG^jr q)F = ieq^F^^{qq) . 



(50) 

Also, in our calculation, we need the interaction part of the quark and gluon background 
fields, 

9sXla{x)xlp{^)[G^,v\cd = {9sqaa{x)[Gf,u]cdqbl3{0))F,(yp , (51) 
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and it leads to the following equation: 



32 

X 



t ~ 

i 1 - ~ 

1 .2, ^ _ (52) 



6 OPE analysis of the correlator 
6.1 Leading order 

Now we calculate the correlation function of the interpolating fields T]n{x) in terms of the 
OPE. First, we carry out the leading-order calculation for the correlator 

n(x) = {QcJp\T{7]r,{x)U0)}\QcJp)F . (53) 

As in Eq. ( 137|1 . this correlator is decomposed into four correlators. We deal with them 
inclusively by using the following notation: 

llkiix) = {^d'\T{jk{x)]im\QcJp)F , {k,l = 1,2) . (54) 

In Figs. [TjlHl the diagrams which contribute to the correlators are illustrated. We denote 
each contribution to the correlators by the upper indices, i.e., Il^^\x) or n^P(x) with 
1 = 1,2, 3. From now on, we use the following abbreviation: 

S'^ix) = CS''^ix)C^ , (55) 

with C the charge conjugation matrix. 

Let us begin with evaluating n^;^(x). Each n[,^'*(x) is expressed in terms of the prop- 
agators Slf^{x) as 

nff(x) = Aeabcea'b'A-Tr[^,Sl,ixh,Si,ix)]^^^^^^ , 
n^^^ix) = 4eabcea'b'A-TTb,SUx)Si,{x)]SU^^^^ , 
n(^(x) = 4e,bcea'b'A-TT[SUx)j,SU^)]j,S^^^^ , 
ng(x) = ^e,bcea't'A-nSL'{x)Si,{x)]j,Si,{x)^^ . 

(56) 
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Figure 1: Diagram which yields the leading contribution without emitting either gluon or 
photon. 



When the propagators include neither photon nor gluon emitting term, the expressions 
reduce to 

U^^^ix) = 24{Tr[755"(a;)755'^(a:)]^'^(x) + 5'^(x)755"(x)755'(x)} , 
u'Six) = 2^TT[^,S-{x)S''{x)]S'{xh, + S'{x)S-{xh,S'{xh,} , 
n«(a;) = 24{Tr[5"(x)75^'(x)]755'^(a:)+755'^(a:)755"(x)5^(x)} , 
ng(x) = 2A{TT[S-{x)S'{x)]^,S''{xh, + ^,S\x)S-{x)S'{xh,} , (57) 

where 

Sl,ix) = 5a,S\x) . (58) 

As discussed in Sec. [3l we focus on the terms proportional to {F ■ cr, ^}. For this reason 
we extract only the terms including {F ■ cr, ^}. They are found to be 



^^^^^^^ ^l^'^^^'^T^^^ ' ^^^}[(^(^dmdpd - eumupu)xG 



uPu Cit mdPd)OQ + {euTTidPu - edmuPd)OG] 
+ {Qdd - du) + (ft - ]:i){&eddd - e„J„)] , 



nS^(a;) =n«(x) 



i 1 ~ 



167r4 (x2)3 

+ {2dd - du) + (ft - ^0(2ed(i<i - Cudu 

nS (a;) (x^^^ ' ^' [{^.edTTidpd - eumupu)xQ 

+ xiieurridpd - edmupu)0Q + {ediriupd - eumdpu)OG} 

+ {Qdd - du) + (ft - ]^C){Qeddd - e«(i„)] . (59) 
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Figure 2: Diagram which yields the leading and the next-to-leading order contributions 
with emitting a photon. Those contributions vanish when P = +1. 



Thus, Il^^\x) is given as 

i 1 ~ 

uPu (^u 

+ 2(3 + 2/3 + 3(3^)dd - (1 + Pfdu 

+ {k- i^) {2(3 + 2(3 + 3(3')eddd - (1 + /3)'e A}] , (60) 

and for /3 = +1, the above expression reduces to 

i 1 ~ 

Il'-^\x) = — {qq)j-^{F ■ a,:jt}[{Aedmdpd - eum^Pu)x0 

+ {Add - du) + {k- ^0(4ed(id - eJu)] . (61) 

Next, we evaluate 11*^^^ Here, we again use the expressions in Eq. ( 1571) . while one of 
the propagators in each correlator is 5'^(x)|i photon in this case. The result is given as 

nf2^(x) =n(?(x) = -^M<ii)j^3i^-^^^} ' 

ng^(^) =^M'ii)j^3i^ ■ ^'^} ' (62) 

and they lead to 

n^'H^) ~, [2(1 - (3)dd - (1 - (3')d.] {m)j^,{F -a^}. (63) 
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Figure 3: Diagram which yields the leading and the next-to-leading order contributions 
with emitting a gluon. Those contributions vanish when /3 = +1. 



Therefore, we find that 11*^^^ vanishes when we take /3 = +1: 

n(2)(x) = . 



(64) 



Finally, we study 11^'^^ In this case, we use the equations in Eq. ( I56l) . By using 
Eq. ( I52l) . we find the resultant expressions as 



nS?(x) 



nS?(x) =n(?(x) 



_ 1 _ _ 1 



[X 



IGtt 

Thus, n(3) (x) is given as 



eddd{K - - {cddu - eudd){K + i^) 



{qq)j^{F-a4} 



X 



[-(1 - [k - ]^Oeddd + (1 - (/€ + (eA - eJd 



Again, the correlator turns out to vanish for /3 = +1: 

n(=^)(x) = . 



(65) 



(66) 



(67) 



Taking the above discussion into account, we conclude that the correlator I{{x) is given 
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as 

n(^) =i^(^"^)(^{^--'^} 

X [(1 + pf{4:ed'mdpd - eumupu)xO 

+ (1 - (5^)x{{(^dmuPu - eumdPd)OQ + {curndPu - edmuPdjoci 
+ (10 + Q(5^)dd - (3 + 2/3 - 

+ d42[(3 + 2/3 + 3/3^) - i(l - /3)2]e4«: - 1^) - (1 - /32)e„(«: + 1^)} 
+ J4(l-/32)e,(«:+^e)-(l + /3)2e4K-ie)}] , (68) 
and its Fourier transform is 



n(.)^./.^..^-^n(.) 



^ :(gg)logf-r^){F-a,^} 



X [(1 + (3Y{Aedmdpd - eumupu)x^ 

+ (1 - /3^)x{(edm„p„ - eumdPdWq + {e^rndpu - edmupd)OG} 
+ (10 + 6/32)drf - (3 + 2/3 - P^)du 

+ dd{2[{3 + 2/3 + 3/3^) - 1(1 - py]ed{K - 1^) - (1 - /3')e.(/^ + 

+ d„{(l-/32)e,(/t + iO-(l + /3)'e4'^-^0}] ■ (69) 

Here, a certain ultraviolet mass scale A is introduced, though it is irrelevant to our final 
result. When one sets /3 = +1, this expression reduces to 

n(g)(°P^) =Y^(gg) log(^) {F ■ i) [i^^dmdPd - eum^pu)xO 

+ {Add - du) + {k- ]^i){Aeddd - eJu)] ■ (70) 

Equation is corresponding to Eqs. (9-12) in Ref. [IS]. After taking /3 = +1, we 
find that the CEDM contribution, i.e., the last term in Eq. fl70p differs from those in the 
reference by a factor of 4. In addition, the sign in front of ^ is opposite to the one in 
Ref. [15]. 
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6.2 Next-to-leading order 

Figures [2] and [3] yield the next-to-leading order (NLO) contributions. By using the prop- 
agator given in Eq. (jUj) , we evaluate the contribution by the diagram in Fig. |2] as 



nl2(a;)NLO = n5i^(x)NLO 



X [eamapdO + {eurriupd - edmdpu)dG + {edUidpd - eumupu)OQ\ , 

(a;)NLO 
i , > 1 



Qo 4(^^)7iT^log(-^/fl^^){^ ■ <y4}(^dmdpdd , 
ng(x)NLO = ^,{qq)7^^og{-K]^x^){F-a4} 



X [edrridPdO - {euniuPd - edmdPu)OG - {ednidPd - eumuPu)OQ] , (71) 



and therefore, n^^^(x)NLO is found to be 



+ (1 - /3^){(e„m„pd - edmdPu)OG + {ednidPd - e«m„p„)6'Q}] . (72) 



n(^n^)NLO =^,{qq)r^^\og{-k]^x^){F ■ a,^}[(l - /3)2e,m,p,^ 



The gluon contribution illustrated in Fig. [3] is also calculated by using the propagator 
displayed in Eq. ( 145|) . The resultant expressions are 

nSf(x)NLO = —^{qq)j^^og{-A^jjix'^){F ■ a,:jt} 

X [edrudpd (ft; +^^)9 + {cdrriupd - eurridPu) (/« - 6*0 

X 

+ {eumdPd - edmuPu){K - ■^^)0q] , 

ni2^(x)NLO = n^i^(x)NLO 

i 1 ~ - 1 

= -^^{^Q)j^^og{-A'jj^x^){F ■ a,:/:}eedmdpd{K + , 

ng(x)NLo = -^(gg)^log(-A?^x^){F-a,^} 

X [edrridPd (^^ + ^0 ^ ^ {edrriuPd - eumdPu) - 6*^ 
-{euTridpd - edrriupu) {1^ - ^O^q] ' 



(73) 
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and then, 

n(=^)(a:)NLo = - ^(gg)^log(-A?^x2){F ■ cr,^}[(l - pyee,m,p,{K + 

+ (1 - f3'^){{edmuPd - eurridpuWc + {eumdPd - edm^pu)6Q}{K - ^Ol • 

(74) 

From the results in Eqs. (172]) and (1711) . it is found that taking /3 = +1 makes the NLO 
contributions vanish, as mentioned before. Thus, we find that the correlator given in 
Eq. f lTOj) is valid up to the next-to-leading order. 



7 QCD sum rules 

In order to derive the QCD sum rules for the present case, we first extract the coefficient 
functions of {F ■ a, ^} from both the phenomenological and the OPE correlators, fl^P^^"^) 
in Eq. and U'^^pe) -^^ (^^^ respectively: 



^(phen)(g2) = 1 



(75) 



^^°''nQ^)-^(^-g)eiog(|!) , (76) 



with = —q^ and 



e = {Aedindpd - e«m„p„)x6' + (4rfrf - d„) + (k - ]^^){Aeddd - eudu) ■ (77) 



In Eq. f[75|) . we neglect the continuum contribution and think of A as a constant, as 
discussed above. The QCD sum rules are obtained by equating the coefficient functions 
after the Borel transformation, i.e., 

^[^(phen)^g2)] = i3 [C^OP'^) (Q^)] , (78) 

where the Borel transformation of the function f{Q'^) is defined as 

with M so-called the Borel mass. Then, we finally derive the sum rules as follows: 

A^d„m„ - AM^ = -e{qq)Mle^ . (80) 

All we have to do is now reduced to determining the Borel mass M and the coupling A„, 
as well as estimating the parameter A. 
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Figure 4: Dependence of our sum rules on Borel mass M^. Range of M is set to be 
0.5 GeV < M < 0.9 GeV. 



To illustrate the dependence of the sum rules on the Borel mass, we plot M^cm^, which 
is included in the right-hand side of Eq. (IHOj) . as a function of the Borel mass squared 
in Fig. H Here, the range of M is set to be 0.5 GeV < M < 0.9 GeV. From the 
figure we find that the Borel mass dependence of our sum rule is moderate in the range 
of 0.6 GeV < M < 0.8 GeV. 

8 Determination of from lattice 

The low-energy constant A„ determines the normalization of the QCD sum rules so that 
the uncertainties are directly linked to the final result. We extract its numerical value from 
the lattice QCD calculation presented in Ref. [23], in which the QCD matrix elements for 
the proton decay rate are evaluated. In fact, they evaluate a similar quantity for proton, 
though the isospin symmetry allows us to interpret it for the present purpose. 

First, we introduce a generic notation for three-quark operators with an arbitrary spin 
structure: 

O^^' {x, t) = eabc [dl{x, t)CTu,{x, t)] T'd,{x, t) , (81) 

with r and F' arbitrary 4x4 matrices. In the current case, the relevant matrices are 
R = Pr = |(1 + 75) or L = Pi = i(l — 75). Now we define parameters ai and 0:2 as 
follows: 

(0|O^^|iV(p,s))=aiPiM„(p,s) , 

(0|(9^^|iV(p, s)) =a2PLUnip, s) . (82) 
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The phase definition is fixed such that ai and a2 are both real and ai < 0. The parity 
transformation of the above equations imphes that 

(0|O^^|iV(p,s))=-aiP^n„(p,s) , 

(0|O^^|iV(p, s)) = - a^PRUnip, s) . (83) 

The interpolating fields ji and j2 are expressed in terms of the operators as 

ji{x) =2(C^^(x) + C^^(x) - O^^(x) - C^^(x)) , (84) 
j2{x) =2(C^^(x) - C^^(x) + O^^(x) - C^^(x)) . (85) 

Thus their matrix elements between the vacuum and one-particle states are given as 

(0|ji|iV(p, s)) =2(ai - a^Wip, s) , (86) 

(0|j2|iV(p, s)) = - 2(«i + a2)un{p, s) , (87) 

and they lead to 

(0|r/„| Ar(p, s)) = 2 [{ai - aa) - /3(«i + «2)] Mn(p, s) . (88) 

From the equation we may relate A„ with the parameters ai and a2'- 

A„(/i) = 2 [(tti - tta) - /3(ai + "2)] , (89) 

with n the renormalization scale. The parameters ai and 02 at ;U = 2 GeV are estimated 
in Ref. [23] as 

«i = -0.0112 ± 0.0012(stat) ± 0.0022(syst) GeV=^ , (90) 
a2 = 0.0120 ± 0.0013(stat) ± 0.0023(syst) GeV^ . (91) 
For /3 = 1, A„(/i = 2 GeV) is given as 

A„ = - 4^2 

= - 0.0480 ± 0.0052(stat) ± 0.0092(3yst) GeV^ . (92) 

Since the QCD parameters used here is evaluated at /i = 1 GeV, we need to translate 
the above value of A„ into that of /i = 1 GeV. The one-loop correction for A„ is 

A„(^ = 1 GeV) = ( ^4^) ( "* A„(^ = 2 GeV) . (93) 

V as{mc) ) V "42 GeV)/ 

which results in a reduction factor of ~ 0.91. As a result, we obtain 

A„ = -0.0436 ± 0.0047(,tat) ± 0.0084(,yst) GeV^ , (94) 

for /3 = +1. 
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Let us compare the value of A.„ obtained here with those used in the previous works. 
In Ref. [H], for example, they exploit the values for A„ evaluated in Ref. (230 by using 
the QCD sum rules. Two Dirac-7 structures, 11 and ^, provide different sum rules. As 
evaluated in Ref. pi], these two sum rules yield relatively small values for A„; the lattice 
QCD value is several times larger than the values evaluated by using the QCD sum rules. 
The author in Ref. [21] also estimates the error for these values. It is about 30 % for the 
sum rules result, while 20 % for the lattice QCD result. The lattice QCD result might 
have a uncertainty in the chiral extrapolation. Since there is no more guiding principle for 
judging which estimation is valid, we exploit the lattice QCD result in Eq. f lM|) because 
this choice leads to rather conservative constraint for CP-violating sources. 



9 Results 

Now we estimate the neutron EDM by using the results obtained above. First of all, we 
rewrite the sum rules in Eq. flSUj) in a simple form: 

Co + cix = fix) , (96) 

where x = and 

/W-^exp(!^). c..^. (97) 



The right-hand side of Eq. fl96|) describes the behavior of the coefficient function obtained 
from the OPE calculation, while the left-hand side represents the phenomenological one. 
The first and second terms in the left-hand side correspond to the double and single pole 
contributions, i.e., the first and second terms in Eq. fl2^ . respectively. Once given a 
Borel mass point x = M^, one may readily pick out cq and ci from the tangent line to 
the function /(x) at the point. Then, they are expressed as the functions of x as 

Co (a;) =^{mlx - x^) exp(— ) , 

1 Tn^ 
ci(x)=-^(2x-m^)exp(^) . (98) 

From these expressions, it is found that the single pole contribution vanishes at a; = m'^/2. 
Since the parameter A is unknown, this choice of x is favorable in order to estimate the 
double pole contribution. Then, at this point the value of Cq is 

Co = 1.8 X 10-2 (for X = ml/2) , (99) 
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Note that the notation used in Ref. [24l is different from ours: 



A„ = 2Ao = TT^Ao . (95) 

Also, notice that there is some difference between the results described in Ref. 122 and the corresponding 
expressions shown in Ref. [2] . 
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Figure 5: Ratio of single and double pole contributions as a function of Borel mass M^. 
Shaded region illustrates that single pole contribution is more than 30% of double pole 
contribution. 



and therefore, the neutron EDM dn is evaluated as 

=z^^0 = 1.2 X 10-^ e . (100) 

Here we take (qq) = -(0.225 GeV)^ . 

The choice of the Borel mass, = "01^/2, is, however, quite arbitrary, and deviation 
from the above result due to the different choice of the Borel mass should be taken into 
account as theoretical uncertainty. In Fig. O we plot the ratio of the single and double 
pole contributions as a function of M^. From this figure we find that the single pole 
contribution rapidly increases when the Borel mass is varied from = m^/2. We here 
assume our sum rules to be valid within the region of the Borel mass in which the single 
pole contribution is less than 30 % of the double pole contribution. This assumption leads 
to 0.36 GeV^ < < 0.50 GeV^, and in this region, dn takes the following values: 

dn = 1.2 X 10"^ e . (101) 

Here, the lower value corresponds to the upper limit of the Borel mass, and vice versa. 

Next, we discuss the uncertainty of the OPE calculation. In this case, the truncation 
of the OPE leads to the uncertainty. Let us estimate it by evaluating the relative size of 
the higher-order contributions. Among them, the four-quark condensates such as (qqqq) 
are expected to yield sizable contributions, since they are free from loop suppression. On 
the assumption that these contributions vanish when one takes the quark masses to be 
zero, we expect that they are suppressed at least a factor of {qq)^ /M"^ ^ 0.1. Therefore, 
the uncertainty of the OPE calculation is estimated to be C(10) %. 



21 



Taking the above discussion into account, we finally evaluate the neutron EDM with 
theoretical uncertainty as follows: 

dn = 1.2 toi ± 0.1 toi X 10"^ e , (102) 

where the first uncertainty stems from the phenomenological calculation while the second 
one comes from the approximation in the OPE. We also include uncertainties originate 
from those in A„ (See Eq. (IM|) .). which is indicated by the third error in the above 
equatioij^. After all, it is found that there is almost C(l) factor of uncertainty in our sum 
rule calculation. 

Let us compare the result with those obtained by using the values of A„ calculated 
with the QCD sum rules. In Ref. |22], the authors adopt = 1.05/(27r)^ GeV^ for 
/3 = —1 from the QCD sum rules derived for the nucleon mass. The A„ for /3 = — 1 is 
equal to that for /3 = 1 in the non-relativistic quark limit. In Ref. [22], it is shown that the 
neutron EDM prediction in the non-relativistic quark model with the SU(6) spin-flavor 
symmetry is derived in the QCD sum rules using the value of A„. By substituting the 
value into Eq. fllOOp . one obtains dn = 3.3 x 10~^ 0. The "realistic" value of A„ evaluated 
by the QCD sum rules in Ref. [21i, A„ ~ 0.022 GeV^, leads to a slightly larger result: 
dn = 4.6 X 10~^ 9. Thus, the overall factors of dn of these results are several times larger 
than that of our result. 

For one's convenience, we substitute the numerical values for the QCD parameters in 
Eq. ([MD. Here we take = 0.8 GeV^ , x = -5.7 ± 0.6 GeV-^ ^ = -0.74 ± 0.2, and 
K = —0.34 ± 0.1 Then, with those parameters the center values, we find 

dn = 4.2 X 10^^^^ [e cm] + 0A7dd - 0.12rf„ + e(-0.18J„ + O.lSda - 0.0084) • (103) 

The contributions from 6 and the quark CEDMs to dn may be changed furthermore by 
about ±10%, mainly due to the theoretical uncertainty of x- 



10 Under the Peccei-Quinn symmetry 

It is known that 0(1) 6 induces too large neutron EDM, the strong CP problem. The 
Peccei-Quinn (PQ) symmetry is one of the solutions for the strong CP problem. If the 
PQ symmetry is introduced, 6 vanishes dynamically. However, if the quark CEDMs are 
non- vanishing, a linear term is induced to the axion potential j27j . 

V{9) = K'e - , (104) 

where K is the topological susceptibility 

K = —i lim i I d'^xe 
k^o J 

^ We approximate the error in A„ as the r.m.s. of the statistical and systematic errors displayed in 
Eq. dMl). 



(105) 
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and K' is calculated by 

K' = -z limz I d'xe'''^T[^GG{x)'- d,qg.s{G ■ a)75g(0)]) • (106) 

q=u,d,s 

Minimizing the axion potential, an effective 6 term is induced from the quark CEDM, 

4. = = !f 5: A . (107) 

K 2 ^-^ rrig 

q=u,d,s ^ 

Taking the induced 6 term into account, the neutron EDM in the presence of the PQ 
symmetry is estimated as 

<Q = 1.2 t'oi ± 0.1 X 10-1 , (108) 

where 

= Md-du + i-^x + f^-^U {^edd - eudu) . (109) 

The contributions from the strange quark CEDM are cancelled in the presence of the PQ 
symmetry [15]. Again, we substitute the numerical values for the QCD parameters as 
presented in the previous section. The result is 

= 0A7dd - 0.12du + 6(0.354 + O.Udu) ■ (HO) 



11 Conclusion and discussion 

We have studied the neutron EDM induced by the CP-violating interactions up to the 
dimension-five operators. In order to derive the relation between the CP-violating inter- 
actions and the neutron EDM, we have used the QCD sum rule technique. There are 
several phenomenological parameters to estimate the relation numerically. Pospelov and 
Ritz also analysed the neutron EDM using the QCD sum rules [TH [T5] and they deter- 
mined the low-energy constant A„ within the framework in the QCD sum rules. On the 
other hand, we have extracted the A„ parameter from lattice calculations. This approach 
allows us to reduce a theoretical uncertainty and leads to a conservative constraint on 
the CP violations. Our result is about 70 % smaller compared with the one obtained by 
Pospelov and Ritz. There still remains a sizable uncertainty resulting from the QCD sum 
rules itself due to a choice of the Borel mass scale. We have estimated the uncertainty 
from the Borel mass scale assuming that the single pole contribution is less than 30 % 
of the double pole contributions. This assumption leads to the theoretical error of about 
0{1). 

Finally, we briefly comment on the contribution of the CP-violating dimension-six 
operators. Among the dimension-six operators, the following operator, called Weinberg 
operator. 
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might be comparable to the quark EDM and CEDM contributions since they suffer from 
the chiral suppression. The other CP-violating dimension- six operators are effective four- 
quark operators of light quarks, which are negligible in the neutron EDM in typical high- 
energy models since the Wilson coefficients are suppressed by the quark masse^. In our 
QCD sum rule calculation it is found that the contribution from the Weinberg operator 
is 0{{qq)'^) and thus sub-dominant. There are a lot of discussions on the significance of 
the Weinberg operator [22], though no consensus has been reached yet. 
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""if the CP-violating four-quark operators include heavy quarks, the CEDMs for the light quarks and 
the Weinberg operators are radiatively induced by integration of the heavy quarks as shown in Ref . [28j . 
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Appendix 

In this Appendix, we list some techniques which we use to carry out the calculation. 



A Quark condensates on the CP- violating background 

In this section, we discuss the effects of the CP-violating interactions on the quark con- 
densates as well as on the quark and gluon background fields. We begin by estimating it 
for the generic quark bi-linear condensate (0|grg|0), with F a 4 x 4 constant matrix for 
which the quark bi-linear qVq is an Hermitian operator. Then, by using the results ob- 
tained there, we derive the relations between the quark condensates and the background 
fields. 

First, we evaluate the quark bi-linear qVq on the CP- violating background. The con- 
tribution of the QCD 6 term at the leading order is evaluated as [30] 

(0|gFg|0),^ = i j d^a;(0|T{g(0)Fg(0)^G^G';:,(a;)G'^'^'^(a;)}|0) + 0{e') . (Ill) 

Substituting Eq. ([5]) into the above expression, we obtain 

mTq\Q)ea=i j <fx{^\T{q{Q)Vq{Q)^-f\d^J,,{x)-2i m,p,q{x)j,q{x)]}\0) + 0(9') , 

q=u,d,s 

(112) 

with pq = Oq/Oq. The first term in the equation is calculated with the aid of the integration 
by parts: 

d'x{Q\T{q{Q)Tq{Qf-fd^J,,{x)m = y"rf^x(0|[J°(x), g(0)Fg(0)]5(x°)|0) 

= ^P.(0|g{75,r}g|0). (113) 

For the second term, we insert the intermediate states and keep only the contributions of 
the one-particle states of the pseudo Nambu-Goldstone bosons vr^ and r^^: 

-^ Yl j d'^{0\T{q{0)^qiO)9GmqPqq{x)i^,q{x)}\0) 

q=u,d,s 
Q 

'ruupu - 'mdPd){qq){0\qTq\7T°) 



Og 



[niuPu + rridPd - '2msPs){qq){0\qTq\r]^) , (114) 



V 

where is the pion decay constanl|f|, and and rrir^ denote the masses of 7r° and ?7°, 



^ We use the PC AC relation for 7r°, 

d^J^ix) - -UmlTTix) , (115) 

and a similar relation for 77". 
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respectiveljl^. As a result, we obtain 



(0|grg|0).o=^P.(0|g{75,r}g|0) 



{muPu - mdPd){qq){0\qTq\7r^) 



0G 

Q 

^rriupu + nidpd - 2msps){qq) {0\qTq\r]°) . (116) 



Other contributions also may be evaluated through a similar procedure. For the 
contribution of the 75-mass terms, 

{0\qTq\0)g = - -^-^(m„p„ - mdpd){qq){0\qTq\7r°) 

{muPu + mdPd-'2msPs){qq){0\qTq\T]^) , (117) 



V 

while for the contribution of the quark CEDM terms, 

1 TTl^ ~ ~ 

(o|grg|o),cEDM =:r^^K - dd){qq){mq\A 

1 Tfl^ ^ ^ ~ 

+ -^{du + dd- 2ds){qq){Q\qTq\r,') . (118) 

Furthermore, it is found that the quark EDMs induce no contribution. 

Taking all of the contributions into account, we obtain the CP- violating contribution 
to the quark condensates as 

(0|gr9|0)gp=^P,(0|g{75,r}g|0) 



1 

+ 



+ 



1 



m5 



2 {du - dd) - 9{muPu - rridPd) 



mo 



— (c?« + dd- 2ds) - 9{rnupu + rUdPd - 2msp. 



{qq){0\qTq\A 

{qq){0\qTq\r]') 



p,(0|g{75,r}g|0) . (119) 



Here, the second equality comes from the conditions in Eq. ( flOl) . Thus, the choice of the 
quark mass phases in Eq. ( |T2l) reduces the contribution of the CP-violating interactions 
to the vacuum condensates into a quite simple expression. 

Note that the CP-violating contribution to the quark condensates vanishes in the 
basis where the 6 term is completely rotated out into the imaginary mass term. Thus, the 



^ The effect of the ir^-rp mixing is suppressed by a small factor of (m„ — m^jms and we ignore it for 
brevity. 
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choice of this basis, which is often adopted in the chiral Lagrangian approach, simphfies 
the calculation. In our paper, however, we remain in a general basis in order to display 
each contribution explicitly. 

Next, we discuss the way of translating the quark and gluon background fields into 
their condensates. To begin with, we consider a single quark line xlai^)Xbi3i^)- ^^i^ 
case, it is related with the quark condensate as follows: 

xLix)xUO) = {^Maa{x)qtm^^)F . (120) 
Using the Fierz identity, the right-hand side of the expression leads to 

6„ 



{^qip\qaa{x)qb0{O)\^qip)F = 



(g(0)g(x))F,QP + 75(g(0)75g(a;))F,QP 

(121) 



+ l(^''"{qiO)cr^,q{x))F,qp 



These quark condensate terms are evaluated by conducting the short-distance expansion 
of the quark field in the Fock-Schwinger gauge as 

g(x) = g(0) + x'^D^g(O) + . . . . (122) 

We note that in this gauge one does not need to care about Wilson-line operators for the 
quark fields. (See Appendix O) 

In the case of CP-even vacuum, the Lorentz and CP invariance of vacuum tell us that 

{qi^q) = {qi^tq) = {qit^i5q) = o . (123) 

On the other hand, on an electromagnetic background, {qu^yq) may have non-zero VEV 
proportional to the electromagnetic field strength F^y. The electromagnetic field depen- 
dence for quark condensates is given as 

(qq), (124) 

where Xq is called the quark condensate magnetic susceptibility [22]. Similar parametriza- 
tion is used for the condensates including the gluon background field: 

gs{qG^yT^q)F = ^qF,y{qq) , (125) 

'2gs{qinGiyT^q)F = i^qF^Am) ■ (126) 
As in Ref. [22], we assume Xqi ^q and to be proportional to the quark charge: 

Xq = egX, f^q = ^gK, = Cg^ . (127) 

This assumption corresponds to neglecting of the closed-loop contribution with gluon 
exchange. 
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Now let us consider the effect of the CP-violating interaction in Eq. ([2]) to the quark 
condensates. By using Eq. f lll9p and the expansion in Eq. f ll22p . we evaluate each quark 
condensate on the CP-violating background as follows (with omitting the subscriptions 
F and QP for simplicity as long as it is not confusing): 



(g(O)g(x)) = {qq) , 
(g(0)75g(a;)) = {.qim)qp = iOcpqiqq) 



(128) 
(129) 



I _ I 



+ -x'^{q[P,a,,]q) 



-mgXf,{qq) 



(130) 



where we use the classical equations of motion in the quark condensates and move the co- 
variant derivatives with help of total derivative. The validity of this procedure is discussed 
in Appendix |Bl Furthermore, 

{Q{^hp.75q{x)) =x''{q-i^D^-f^q) 

1 _ i _ 

= ^x''g^u{qJpl5q) + ■^x''{q[p,cXf,^]-f5q) 

t — ~ % 1 ~ 

= 2'^q(^qX{0pqF^,y + Ff,^)x'' (qq) + -{dq -C]eqdq)x''F^^{qq) , (131) 



and 

{Q{0)(^f^uq{x)) ={q(^^uq) = {q(^t,uq) CP even QP 
=eqx[Ff,^ - OGpqFf,u]{qq) ■ 
Taking the above discussion into account and using the relation, 



F^,x^'Y'y5 = +^{F-a,:/:} , 



and 



F ■ a = iF • (J75 , 
we finally obtain the expression for the single quark line as follows: 

xlMxlp{^) = - ^ (1 + iOGpql5)o,p {qq) + j^Sab^^a^mgiqq) 



(132) 

(133) 
(134) 



12 
96 
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OrnqPgCqX + dg+{K,- -O^qdiq 



{F ■ a,:jt}a(3{qq) 



{qq) ■ 



(135) 
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Lastly, we evaluate the interaction part of the quark and gluon background fields, 

9sXLix)xlp{^)[Gt^u]cd = {gsqaa{x)[Gf,^]cdqbi3{0))F,qp ■ (136) 

Again, we use the Fierz identity and the short-distance expansion of the quark field, and 
through a similar calculation, we find the following results: 

9sXLix)xl0{O)[Gfiu]cd = - ^{SadSbc - ^SabScd){qq) 

i 1 - ~ 



12 



-rrir 



(137) 
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B Equations of motion 

Let us discuss the validity of using classical equations of motion for quark condensates. 
We investigate the following quantity: 



(0|gr(^^-m,)g|0), 



(138) 



where F is a (c-number) 4x4 matrix. The subscript indicates that this quantity is 
evaluated in the 9 vacuum. One may readily generalize the discussion here for the case 
with other CP-violating sources. The discussion presented in this section is based on 
Ref. [31] • 

First, we define the generating functional Z[ri\Q on the same background: 



VqVq exp 



Here, the Lagrangian density C is 



C = qiilp — mg)q . 



(139) 



(140) 



The functional derivative of the generating function with respect to the function tj yields 
Eq. (fT38D . that is, 

{O\qr{tp-mq)q\0)eoc^^^'^^' 



571(0) 



(141) 



77=0 



Now we replace the integration variable q with a new integration variable q' as 

q ^ q = q — rjqT . (142) 
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Since this step does not change the integral, then we obtain 



VqVq 



Det 



5^ 



exp 



d^xC + 0{rf 



(143) 



where the inverse of the Jacobian comes from the transformation of the measure for the 
fermionic variable. 

Next, we evaluate the Jacobian in the expression above. Since we are interested in the 
first order derivative of the generating function, we expand the Jacobian in 77 and keep 
only terms linear in 77. 



6qa{x) 



S^pS^x -y)~ v{y)T^p5\x -y) + Oiv') 



Using the identity 



DetM = exp Tr In M 



we readily obtain the Jacobian as 



6q 



exp 



Tr(r) / d'^xri{x)5'^{x — x) 



(144) 
(145) 

(146) 



It is found that if the trace of the matrix F is non-zero, the Jacobian yields a singular 
factor, while if it vanishes we need careful treatment for evaluating this term. So, in the 
following discussion, we divide F into two types; one is the term proportional to the unit 
matrix and the other is the traceless part. 

First, we consider the case F oc 1. Using Eqs. fll4ip . (11431) . and (I146p . we obtain the 
following equation: 

{0\qT{tp - mg)q\0)e = -2Tr(F)5^(0) . (147) 

Once you carry out the normal ordering for the composite operator qT{ip — mq)q, the 
singular factor in the right-hand side vanishes. Thus we conclude that after normal 
ordering. 



(0|gF(z|)-m,)g|0)e = 
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when F oc 1. This equation implies that we may use the equations of motion for quark 
condensates in this case. 

Next, we shall turn to the traceless part. In this case the Jacobian in Eq. f ll46p is 
written in terms of the anomaly function defined as 



A{x) = 2TT{r)6\x -x). 
With this function, Eq. f ll46p leads to 

Det"^ 



Sq^ 



exp 



d xri{x)A{x) 



(149) 



(150) 
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The usual analysis for the chiral anomaly tells us that the function A{x) in Eq. (11491) does 
not vanish only for the case F = 75. Thus, if F 7^ 75, Eq. (11481) is satisfied. When F = 75, 
on the other hand, the function A{x) is evaluated as 

Ax) = ^G%G^>^^ . (151) 

From Eqs. (fTiTj) . (fTi3|) . and (1T50|) we eventually find that 

(0|g75(^^ - m,)g|0), = -^(0|G^,G^'^10), . (152) 



57r 



(155) 



This expression is simplified via the axial current anomaly equation: 

5^( Wsg) = 2zm,g75g + "^G^^G^^^ . (153) 

Then, Eq. (1T52|) leads to 

{0\q^MO)e = ^(0|9'^(g7M75g)|0)e = . (154) 

Therefore, we are not able to use the classical equations of motion for quark condensates 
in this case. 

As a result, we find that 

{0\qTi]pq\0)e 
^ f {0\qTmgq\0)g (for F = 1, 7^, 7^75, a 
~\0 (forF = 75) 

Also, its conjugate leads to 

- {0\qt^Tq\0)e 

{0\qmgTq\0)f) (for F = 1, 7^, 7^75, cr^^,) 
(forF = 75) 

Before concluding the section, we add a comment on the condensate of the total 
derivative terms. As we have already conducted in Eq. (I154p . the Lorentz invariance of 
vacuum implies that condensates of the divergence of quark bi-linear always vanish, i.e., 

{0\d>'{qT,q)\0) = 9'^(0|(gF^g)|0) = , (157) 

with F^ a constant matrix which transforms as a vector under the Lorentz transformation, 
such as 7^, 7/^75, or so. On the other hand, the total derivative of the quark bi-linear is 
written as 

d'^iqT^q) = (9^g)F^g + qT^{d^q) = qt'^T^q + qV^D^q . (158) 

Thus we find 

(0|g^'^F^g|0) = -{qT^D^q\Q) . (159) 
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C Wilson Line in Fock-Schwinger gauge 



Quark fields q{x) are always accompanied by an appropriate Wilson-line operator in order 
to compensate the different gauge transformation property of the quark fields at different 
space-time points. In the Fock-Schwinger gauge, however, one may always choose a 
particular path which makes the Wilson-line operator equal to identity [32]. We show 
this statement in the following. The Wilson line is written as 



Up{x,0) = P |exp 



jexp Ws 



\s^^^G^{x'{s))T^ 
ds ^ 



where 



x'(0)=0, x{l)=x 



(160) 



(161) 



and P denotes path-ordering. This operator depends on the choice of the integration 
path. Here we take a path such that 



x'{s) = sx 



(162) 



Then, 



Up{x,0) = pjexp 



Ws 



dsx>'G^{sx)T' 



(163) 



In the Fock-Schwinger field, the gluon field is expanded as 



G„{x) = ATa;^G,^(0) + ;^a;°a;'^(D„G',^(0))+^a;"a:V(D„D;,G',,,(0)) + --- • (164) 



2 • 0! 



3-1! 



4-2! 



Inserting this expression into Eq. f ll63p . we readily find that all terms in the exponential 
vanish due to the antisymmetric property of the gluon field strength tensor. Therefore, 



Up{x,Q) 



(165) 
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